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We investigate the vector-current correlation (VCC) in the presence of a strong external 
magnetic field (B = Bqz) at low temperature (T < T^) with C- and CP- violations, indicated 
by the nonzero chiral-chemical potential {^ix 7^ 0)i the chiral-magnetic effect (CME). For this 
purpose, we employ the instanton-vacuum configuration at finite T with nonzero topological charge 
{Qt 7^ 0). From the numerical calculations, it turns out that the longitudinal component of the 
connected VCC is liner in Bo and shows a bump, representing a corresponding vector meson at 
IQI — (300 ~ 400) MeV for T — 0. The bump becomes enhanced as T increases and the bump 
position shifts to a lager \Q\ value. In the limit of \Q\ — >■ 0, the transverse component of the 
connected VCC disappears, whereas the longitudinal one remains finite and gets insensitive to Bo 
with respect to T, due to diluting instanton contributions. Considering a simple collision geometry 
of HIC and some assumptions on the induced magnetic field and screening effect, we can estimate the 
charge separation (ChS) as a function of centrality using the present results for VCC. The numerical 
results show a qualitative agreement with experiments for the Au+Au and Cu+Cu collisions. These 
results are almost independent on the source of CME, instanton or sphaleron, as long as the CME 
current is linear in Bo. 
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I. INTRODUCTION 

QCD at finite T is one of the most important and intriguing subjects for decades. Along with the energetic progresses 
achieved in heavy-ion coUision (HIC) experimental facihties such as the relativistic heavy-ion coUision (RHIC) at BNL, 
theoretical studies to understand QCD at finite T become far more exciting subjects than ever before, especially for 
the vicinity /i ~ and T ^ 0, resembling the early universe. Beside the fact that the bulk properties of HIC can be 
interpreted well by relativistic hydrodynamics [Tl[2], the microscopic features still seems to be governed by QCD. 

Recently, it was proposed that P- and CP-violations due to the nontrivial QCD vacuum can be measured experi- 
mentally by seeing event-by-event charge separations (ChS) in HIC experiments [3HSj- Briefly, this phenomenon can 
be understood as follows: 1) A non-central collision of two heavy nuclei can generate a strong magnetic field, which is 
perpendicular to the reaction plane. 2) If there is nonzero chiral-chemical potential fi^ ~ Nn — 7^ proportional 
to nontrivial topological charge Qt, indicating tunneling between QCD vacua and P- and CP- violations as well, 
3) quarks with different chiralities (or helicity for massless quarks) move in different directions along the magnetic 
field. 4) Simultaneously, according to the nonzero /i^, there appear electric currents, produced by left-handed or 
right-handed quarks, resulting in ChS in the measured particles. This is usually called the chiral-magnetic effect 
(CME). 

Interestingly enough, it was reported that this intresting phenomena, CME was indeed observed by the recent RHIC 
experiments by STAR collaboration [71 15] . Since ChS is a consequence of event-by-envet fluctuations, if one average it 
over a space-time volume, it disappears. Thus, a non-vanishing correlation between two particles with electric charges 
a and b, measured in different azimuthal angles (Ac/)), was suggested [HIU]: 

Charge separation (ChS) = ((cos(A0a -t- A0f,))), (1) 

where ((•••)) denotes an average over the space-time volume. According to this, the experiments observed sizable 
strengths for ChS for the Au-|-Au and Cu-|-Cu collisions in non-central collisions. Theoretical studies also have been 
done energetically for this direction in various approaches: model- independent analyses on CME jTDHT^ . lattice QCD 
simulations [T31 [Tl] , effective QCD-like models , and gauge-gravity duality inspired models [TOH^ , and so on. 

In what follows, we briefly introduce these theoretical works. In Refs. [10], it was shown that the induced current due 
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to CME along the magnetic field can be derived as: 

(J||)B..,=iV,7V,g^, (2) 

where the subscript \\ stands for that the current is parallel to the magnetic field B = Bqz. Here, we set the quark 
electric charge unity for brevity. As indicated in Ref. [lOj and references therein, this expression is a very general 
consequence of the axial anomaly in QED. The lattice QCD simulations indicated that the longitudinal component of 
the current is much larger in comparison to the transverse one. In addition, it was shown that CME becomes insensitive 
as T increases. These interesting results were also supported and reproduced by the instanton-vacuum configuration 
model at low T [151 US]- On top of the (partial) chiral-symmetry restoration, deconfinement phase of QCD was 
also studied with CME employing the Polyakov-loop-augmented Nambu- Jona-Lasinio (PN JL) model [T71 [TH] . The 
chiral-magnetic spiral (CMS) was also investigated with the angular momentum produced in the non-central HIC |22| . 
Some controversial consequences were reported in Refs. [TM^ . in which the gauge-gravity dual model (AdS/CFT or 
AdS/QCD based) was exploited. Especially, Ref. [20] suggested the negligible CME in the strong-coupling limit. An 
interpretation were given for this dispute with discussions on back-reactions in Ref. [53j. Although the experimental 
observations and theoretical estimations seem promising, it is an urgent task to confirm whether the experimental 
results are genuine consequences of the nontrivial QCD vacuum effect, i.e. nontrivial correlations between two different 
vacua Qt 7^ by the possible sources for CME: instanton for low T [MH2S] or sphaleron for high T [?7l 128] . 

Considering all of these recent progresses for CME, in the present work, we want to study the vector-current 
correlation (VCC) in the presence of an external magnetic field and nonzero Qt at low T, T < T^, which is the 
critical T for chiral restoration, employing the instanton-vacuum model as in the previous work [151 . This model 
is characterized by two phenomenological parameters, i.e. average inter- (anti) instanton distance _R « 1 fm and average 
instanton size p ~ 0.3 fm |25j . We take into account T modifications on these quantities by introducing the Harrington- 
Shepard caloron, which is a periodic instanton in the Euclidean temporal direction with trivial holonomy, but without 
confinement [291I30| . As a consequence, we obtain T- and momentum-dependent constituent-quark mass, which plays 
a role of a natural UV regulator. Since the vacuum-expectation value (VEV) of the induced CME current, {J\\,±} B,fi^ 
was already computed in the previous work [151 116j . we focus to compute the connected VCC in the present work. 
Using the previous and present results for VCC, we also want to estimate ChS by the following relation [5l|9|: 

ChS EE ((cos(A0, + A0b))) « (J||)'b,^, - {J±fB.^, +X1|.B,M, - X±,B.,^., (3) 

" V ' ^ V ' 

Disconnected Connected 

where the first and second two terms in the right-hand side come from disconnected and connected quark-loop 
diagrams of VCC, respectively. Hence, the disconnected VCC can be represented by the squared VEV of the induced 
current, whereas the connected one by the susceptibility x ^ d^S/dA^ . The disconnected diagram contributes to 
CME primarily, while the connected one corresponds to backgrounds [9]. It is worth mentioning that ChS can be 
sorted to two different kinds: (a, &) — (±,±) and (a, 6) — (±,=f), where ± stands for the electric charges of the 
detected particles at different azimuthal angles At/)^ and A0b in HIC. Considering the geometrical symmetry of the 
non-central HIC, ChS for (4-, -I-) and (— , — ) will give the same results, whereas that for (±, =f) is different from others 
in principle. Thus, one can define two different correlations for ChS as same-charge (SCC) and opposite-charge (OCC) 
correlations, i.e. ((cos(A(/>-|- + /S.(f)±))) and {{cos{/S.(j)± -I- A0;p))), respectively. 

From the numerical calculations, it turns out that the longitudinal component of the connected VCC is a liner 
function in Bq and shows a bump, representing a corresponding vector meson at \Q\ = (300 ~ 400) MeV for T — Q. 
The bump becomes enhanced as T increases and the bump position shifts to a lager \Q\ value. In the limit of \Q\ — 0, 
the transverse component of the connected VCC disappears, whereas the longitudinal one remains finite and gets 
insensitive to Bq with respect to T, due to diluting instanton contribution. Considering a simple collision geometry 
of HIC and some assumptions on the induced magnetic field and screening effect of quark-gluon plasma (QGP), we 
can estimate ChS as a function of the centrality. The numerical results show a qualitatively good agreement with 
experiments for the Au-|-Au and Cu-|-Cu collisions. We observe that the strength of ChS for the Cu-|-Cu collision 
is generically larger than that for the Au-|-Au one as a function of centrality. This tendency is caused by that the 
probability for the domain with Qt created in HIC does not depends on centrality, but is proportional to the 
size of the nucleus inversely [5l[7l|8]. Moreover, the absolute value of ChS for OCC becomes smaller than that for 
SCC, due to the screening effect. We note that these results are almost independent from what the source of CME 
is, instanton or sphaleron, as long as the CME current being linear in iiy,BQ. 

We organize the present work as follows: In Section II, we briefly introduce the instanton-vacuum model and define 
a quantity 5, which indicates the strength of P- and CP- violations with nonzero Qt. We obtain an effective quark- 
instanton action as a function of 6. We define connected VCC and corresponding scalar VCC in the presence of the 
external magnetic field and d in Section HI. In Section IV, T modifications on the instanton parameters R and p are 
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taken into account in terms of the Harrington-Shepard caloron, and the constituent-quark mass is then defined as a 
function of momentum and T. Considering all the ingredients discussed in the previous Sections, we write expressions 
for VCC relating to CME as functions of T, 5, and Bq in Section V. In Section VI, we discuss a simple collision 
geometry of HIC. Numerical results for VCC and ChS are given in Section VII with discussions. The final Section is 
devoted for summary and conclusion of the present work. 



II. EFFECTIVE ACTION FROM THE INSTANTON VACUUM WITH P AND CP VIOLATIONS 

In this Section, we briefly introduce a P- and CP- violating effective action ScS: derived by Diakonov et al. from 
the instanton- vacuum configuration in the large limit at zero T [31 . Employing a dilute grand canonical ensemble 
of the (anti)instantons with finite instanton-number fluctuations, A = — iV_ ^ 0, which corresponds to a CP- 
violating vacuum, but a fixed total number of the pseudo-particles + = N, Scs can be written in momentum 
space with the Euclidean metric as follows: 



?cff 



C 



N+ 
V 



_/V_ uiN 
lnA+ + — lnA_--^(A+ + A_) 

>-|[A+(l + 75) + A_(l-75)]p2(fc) 



^ — im 



(4) 



where we have used Nf = 1 for simplicity. However, the extension to an arbitrary Nf is just straightforward. C 
stands for an irrelevant constant for further investigations, whereas N±/V for the (anti)instanton packing fraction 
proportional to the inverse of the average inter-(anti)instanton distance w (200 MeV)*. A± denotes a Lagrangian 
multiplier, which was employed to exponentiate the 2Nf-\ Hooft interaction in the effective action ^5]. The average 
instanton size in the dilute instanton ensemble is assigned as 1/p ~ 600 MeV, while m indicates a small but finite 
current-quark mass for the SU(2) light-flavor sector (to — > 0). F{k) denotes the quark form factor originating from 
the non-local quark-instanton interactions and is deflned as 



Io{t)Kiit) - h{t)Ko{t) - -^h{t)Ki{t) 



t =^ 



\k\p 
2 ' 



(5) 



where /„ and Kn stand for the modifled Besscl functions. We, however, will employ a parameterization of this from 
factor for convenience in the numerical calculations in what follows. 

From the effective action, we can obtain the following two self-consistent (saddle-point) equations with respect to 
A±: 



95, 



off' 



V 
N± 
V 



\j.mNc 



N, 



(1 ± S)MomN, 



d^k 

W) 

-Nr. 



jTr^ 



(l±75)^^^(fc) 



-|[A+(l+75) + A-(l-75)]i^2(fc) 



Tr 

(27r)4 



i(l±75)(l + <575)'M' 



fc2 -H (1 575)2Af2 



0, 



(6) 



where A± is approximated as Mo(l ± 5) in the last line of Eq. ^ accounting for A <C iV in the thermodynamic 
limit [31]. The momentum-dependent constituent-quark mass is defined as Mk = M^F'^ik) [25]. By adding and 
subtracting the instanton (4-) and anti-instanton (— ) contributions in Eq. (|6]), we arrive at 



N 
V 



mMoNc 
27r2p2 



4iV, 



d*k {1 + 5^)AP 
(27r)4 k^ + {l + 6^)AP'' 



(7) 



A SmMpN, ^ f d^k SAP 

V 27r2p2 ^ V (27r)4fc2 + (1 + 52)^/2- y) 

Taking into account 5 ^ 1, A <C iV, and using Eq. ([s]), we can obtain an expression for (5 as a function of relevant 
parameters: 



, ■ (9) 



FIG. 1: Connected vector-current correlation (VCC) in the leading contribution. The wavy and solid lines indicate the vector- 
particle and quark lines, respectively. 



This equation tells us that S contains the information on the instanton-number fluctuation A at a certain scale p, 
which is about 600 MeV in the present framework for vacuum. Taking into account all the ingredients discussed so 
far, finally, we can write the relevant effective action with A for further investigations: 



ScS 



W) 



rTr, 



cj,7 



In 



K 



(10) 



where the Trc,/,-y denotes the trace over color, flavor and Lorentz indices. 

Now, we are in a position to discuss the relation between the nontrivial topological charge Qt, as a source of CME, 
and the instanton number fluctuation A. According to the axial Ward-Takahashi identity [5], Qt is proportional 
to the number difference between the chirally left- and right-handed quarks, Qt oc Nji — N^. Hence, nonzero Qt 
indicates the chirality flip. Note that, similarly, if a chirally left-handed quark is scattered from an instanton to an 
anti-instanton, the quark helicity is flipped to the right-handed one, and vice versa. This means that the nonzero 
A results in Nj^ — Nj^ ^ 0. In this way, Qt can be considered to be proportional to A: Qt ^ ^ [Ml US]. As a 
consequence, we can study CME using the effective action in Eq. (10) as a function of Qt, more explicitly 5 oc A. 



Since 5 stands for the symmetry breaking, its order of strength must be very small, corresponding to the order of CP 
violation. In Table |l] we summarize the numerical values for the relevant inputs. 



III. VECTOR-CURRENT CORRELATION IN THE PRESENCE OF Bo AND 5 



In this Section, we discuss and derive VCC from the effective action in Eq. ( 10 1. First, we write a local vector-current 
operator in Euclidean space as follows: 

V';,{x)^-zq\x)^J-^q{x), (11) 

where q{x) and r° stand for a quark field and SU(2) isospin matrix, respectively. Using this definition, we can write 
VCC (V°), represented by VEV of the time ordering of the two vector-current operators: 

KM = -J d^xe^'^-^0\T[V^,{x),Vt(0)] |0), (12) 

where superscript a and momentum q denote an isospin index and the momentum of a corresponding vector meson, 
respectively. Taking into account vector-current conservation, one obtains an expression for a scalar VCC from 
Eq. ( 12 ) as follows: 

n;it('?)= k.-^')n'^'('z)- (13) 



R 


P 




Mo 


m 


1 fm 


1/3 fm 


197 MeV 


325 MeV 


5 MeV 



TABLE I: Instanton parameters, constituent and current quark masses in vacuum. 



5 



According to this, one is lead to a simple relation between the tensor and scalar VCC: 



Substituting Eq. (14) into Eq. (12 1, we have the following equation: 

d^xe^''-^0\T[V^,{x)XiO)] |0), 



IT 



(14) 



(15) 



In order to compute VCC in Eqs. (13 1 and (15), we rewrite the effective action with an external-source vector field 



V°,, using Eq. (10): 



^off — — 



(2^ 



Tr,j,^ In - z(l + 5-i^)Mk + T^c.] 



(16) 



where we have set the current-quark mass zero in the effective action, since it is irrelevant in further discussions. K^^ 
indicates covariant quark momentum, gauged by a photon field as fc^ -|- A^, in which we set the quark electric charge 
unity for convenience. The electromagnetic (EM) field configuration is chosen to make the external magnetic field 
along the z direction, B — Bqz, which is set to be perpendicular to the collision plane of HIC: 



A,. 



A 



clas 



A 



flue 



^X2,^a;i,0,0) +(AA4). 



(17) 



Here the fluctuations are assumed to be static and far smaller than the model scale: {\A\, A4) <C A w 1/p. According 
to this configuration, only F12 and its dual F34 will survive in the EM field-strength tensor. Then VCC can be easily 
evaluated by performing functional derivatives of Eq. (16) with respect to V^{x) and V^(0), resulting in 



(o|r[v^(x),v,^(o)] |o) = z 



(18) 



In order to induce an external electromagnetic field to the quark- instanton system as Eq. (17), we have used the 



linear Schwinger method [35]. F°i^(p) for connected \CC is then given by 



i5- 



ah 



(2^ 



rTr 



1 



-7m 



(19) 



Here, we have used Tr[T''r''] — 2(5°''. The gauged quark propagator inside the square bracket of the right-hand side 
can be expanded in terms of (5 and A^, and we obtain 



1 



1^ + /I + i(l + 575) {M + \M{o ■ F)\ 



-i{l + 5-i5)M{K) 



1 



fc2 + (l+J2)7\,f2 

M{a ■ F) + iM{k)-f^K,F^^ - 2iSM^5f< 
fc2 + (l + ,52)M2 



(20) 



where a ■ F = Uf^uF^'^ with the antisymmetric tensor ct^i, = ii'j^jnu — 7!^7/i)/2 and the EM field strength tensor 
piiu _ dfj^Ai, — diyA^. The relevant mass functions in Eq. (20 1 are also defined as: 

Mfc = Mo 



Mfc = - 



(21) 



Using Eq. (20) and collecting the terms of 0{5), the leading contributions for the transverse (/i, v = 1, 2), longitudinal 
(/.t, V = i), and temporal (/i, v = A) components of Eq. ( 19 1 can be obtained as follows: 



^1'U2(P) 
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4:iS''''6BoNcNf 



d^k 

(27r)4 (fc2 + M2)[(A:-p)2+Af2_^] 
d-^k 



(27r)4 (fc2 + M2)[(fc-p)2+M2_^] 



MkMk 



k^ + Ml ik-p)^ + Ml 

MkMk Mk-pMk-p 
k-^ + Ml ^ {k-pY + Ml_ 



ei2-Pfc,(fc-p)p, 
[ei2-PA:,(A:-p)p + fc4P3)] 
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-^44^) 



AiS^^SBnNrN 



f 



(27r)4 {k^ + Ml)[{k - py 



ML,] 



MkMk 



Mk-pM, 



pIVlk-p 



fc2 



Mi 



{k-pY + Ml_^ 



[ei2-PA:,(fc-p)p-fc3P4)],(22) 



where we have used the fact that, due to our choice of the EM-field configuration in Eq. (17), only F12 = Bq survives 



among the field-strength tensor. It is clear that the sign difference in the square bracket in the right-hand side makes 
the longitudinal and temporal contributions -^33 44 are significantly larger than that for the transverse ones, which is 
a general consequence of CME as indicated in Refs. [TS]. The term with e tensor can be evaluated further as 

e^^'^''fco.(fc - p)p = A;3(fc4 - Pi) - k^iks ~ pz) = A;4P3 - ^3^4- (23) 



Finally, we arrive at the following expressions for connected VCC and corresponding scalar VCC using Eq. ( 13 ) as 
follows: 



rf4 



(27r) 



d-^p 
(2^ 



-i(p—q)-x ' 



^;tip) = E-^nt(9) - -\ [2J-4t(<Z) +-^33(9) 



--^44(9)] 



(24) 



n=l 



Since p and k are the integral variables as in Eqs. (22 1 and (24), the interchange k ^ p does not make any difference 

However, as will be 



in§3'(g)l 



inJtCz)!- 



in the absolute values of the longitudinal and temporal components, i.e. |±±33\^yy| — |-i-i44\ 
shown in Section V, they become differing from each other at finite T, due to the breakdown of Lorentz invariance, 
since the fourth component of four momenta becomes periodic in the temporal direction. 



IV. INSTANTON EFFECTS AT FINITE TEMPERATURE 



To investigate the physical quantities in hand at low but finite T {T < T^), we want to discuss briefly how to 
modify the instanton variables, p and R at finite T. We will follow our previous work [33] and Refs. [29l [30] to this 
end. Usually, there are two different instanton configurations at finite T, being periodic in Euclidean time, with trivial 
and nontrivial holonomies. They are called the Harrington-Shepard [3D] and Kraan-Baal-Lee-Lu calorons [5^ I35j. 
respectively. The nontrivial holonomy can be identified as the Polyakov line as an order parameter for the confinment- 
deconfinement transition of QCD. However, since we are not interested in the transition in this work as long as we are 
at relatively low T, we choose the Harrington-Shepard caloron for the modifications at finite T. Note that here are 
two caveats: 1) These modifications are done for a pure glue system without dynamical quarks. Hence, the instanton 
variables may change, if one takes into account dynamical-quark contributions in the instanton distribution function. 
2) Moreover, we assume a CP-invariant vacuum for the modifications, i.e. A = 0, whereas we are interested in the 
physical quantities for A 7^ for CME. Correcting these inconsistencies may give rise to changes in the final results, 
although they seem small considering that the order of CP violation in reality is very tiny, but it must be beyond our 
scope in the present work. Keeping this issue in mind, we write the instanton distribution function at finite T with 
the Harrington-Shepard caloron as follows: 

d{p, T) = Cn^ A^s /3^= ^"-5 exp [-{A^T^ + hnp')p'] • (25) 

■■^ ^ ' 

c 

Here, the abbreviated notations are also given as: 

4 gQ g-1.68aRs iVc 

/3 = -61n[ARsPcut], ;3 = -61n[AR,s(i?)], Cn^ ^ ^^.^ _ (26) 



An, = 



11 



-Nr-l 



7 



21 



b = 



llNr - 2N 



f 



N 
V' 



(27) 



Note that we defined the one-loop inverse charge /3 and /3 at a certain phenomenological cutoff value p^ut £^nd 
(i?) « R. As will be shown, only /? is relevant in the following discussions and will be fixed self-consistently within 
the present framework. Ars stands for a scale, depending on a renormalization scheme, whereas V3 stands for the 
three-dimensional volume. Using the instanton distribution function in Eq. (25), we can compute the average value 
of the instanton size, straightforwardly as follows [24 : 



fdpp^d{p,T) 
Jdpd{p,T) 



AnT^ 



2/37^ 



(28) 
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where ly = {b — 4)/2. Substituting Eq. (281 into Eq. (25), the distribution function can be evaluated further as: 

'1 



d(p,T) =C A^exp [~M{T)p^] , MiT) - ^A^T^ 



(29) 



The instanton-number density n can be computed self-consistently as a function of T, using the fohowing equation: 

n^M{T)^[CT{u)Y , (30) 

where we have replaced NT/V3 — > n , an d r( i^) indicates a T fucntion with an argument ly. Note that C and f3 can 
be determined easily using Eqs. (28) and (30), incorporating the vacuum values of the n and p: C k, 9.81 x 10"'' and 



/3 « 9.19. At the same time, using these results, we can obtain the average instanton size p as a function of T with 
Eq. (|28l). 



Finally, in order to estimate the T dependence of the constituent-quark mass Mq, it is necessary to consider the 
normalized distribution function, defined as follows: 



dNip,T) 



d{p,T) _ pfe-5>i'-(T)exp [-M{T)p^] 
Jdpdip,T) ~ W) 



(31) 



Now, we want to employ the large- A^c limit to simplify the expression of dN{p,T). Since the parameter b is in the 
order of 0{Nc) as shown in Eq. (26), it becomes infinity as Nc — )■ 00, and the same is true for v. In this limit, as 
understood from Eq. (31), (ijv(p, Tjcan be approximated as a 5 function |31) : 



lim dNip,T) ^ 5[p- p{T)]. 
Considering the constituent-quark mass can be represented by |31j 

Mo'xV^ f dppH[p - p{T)] = ^/^)p\T), 



we can modify Mq as a function of T as follows: 



Mn 



Mo{T) 



(32) 



(33) 



(34) 



where we will use AIq sa 325 MeV as done for zero T. The numerical results for the normalized p/po and n/no as 
functions of T are in the left panel of Fig. [2j As shown there, these quantities are decreasing with respect to T as 
expected: decreasing instanton effect. However, even beyond « Aqcd ~ 200 MeV, the instanton contribution 
remains finite. In the right panel of figure, we draw the quark mass as a function of T and absolute value of three 
momentum of a quark |fe|: 



M(|fc|,r) = Mo(T) 



2 + pHT)\k\ 



(35) 



Note that we have ignored the Euclidean-time component of the four momentum by setting ^4 = 0. This tricky 
treatment simplifies the calculations in hand to a large extent, and we also verified that only a small deviation 
appears in comparison to full calculations. Moreover, p in Eq. (35) is now a function of T as demonstrated by 



Eqs. (28) and (34) previously. As shown in the figure, M{\k\,T) is a smoothly decreasing function of T and |fc|, 
indicating that the effect of the instanton is diminished. For more details, one can refer to the previous work [33] . 



VECTOR-CURRENT CORRELATION AT FINITE TEMPERATURE 



In this Section, we briefly discuss how connected VCC in Eq. (24) are converted as a function of T. For this purpose. 



we make use of the fermionic Matsubara formula (Euclidean-time formula) as follows: 



d^k 

J2ny 



F{k) T 



7l — ~00 



d^k 



F{k,Wn), 



(36) 
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FIG. 2; (Color online) Normalized p/po and n/no as a function of T for A'^c = 3 (left). AI as a function of T and absolute value 
of the momentum |fc| (right). 



where the fermionic Matsubara frequency reads w„ = {2n + 1)ttT. Using this formula, the relevant four vectors in 
Eq. m can be rewritten by 



e = 



{k - qf = wl + E, 



(37) 



where k — (0, ksin 6', kcos 9,Wn) with k = |fc| and g = (0, 0, 0, IQj) in the space-like region for the corresponding vector 
meson at rest for simplicity. The energies of the quarks involved are also defined by 



E' 



k^ + {l + d^)Ml El 



2«;„|Q| + |QP + (1 + <5')A42 



2wo\Q\ + \Q\^ + {l + 5^)Ml 



(38) 



Note that we have kept the term 2u;„|(5| for Ei, and assumed that the lowest Matsubara frequency = ttT is 
dominant. Later, we will see that to keep this term is important to have correct \Q\ behaviors of VCC. Similarly, the 
T-and momentum-dependent quark masses in Eq. (21) are also redefined as follows: 



Ma = Mo 



k2p2 



Mb = Mo 



2 + (k2-27«o|Q| + |Q|2)p2 



(39) 



The T dependence of these masses are included in A/q as shown in the previous Section. As understood in Eq. (21 1, 
Ma,b and are also similarly redefined. We write a useful summation identity which will be used to evaluate the 
denominator in the expressions for VCC: 



E 



T 



1 



8TXy3(X2-y2): 



-sech 



Y 
2f 



XY^ - YX^ + 4Tr3cosh^ 



Y 
2f 



tanh 



2T 



TX{X^ -3Y^)smh 



= M{X,Y). 



(40) 



Considering all the ingredients discussed so far, finally, we arrive at expressions for the imaginary parts of connected 
VCC as a function of T, strength of the external magnetic field Bq, and strength of P and CP violations 6: 



Imnii^22('7) 
Imn33(q) 
Imn44(q) 



-2\Q\SBoN,Nf 
-2\Q\6BoNcNf 
-A\Q\SBoN,Nf 



k^dk cos 0d cos 6* 

47^2 

k^dk cos 0d cos 6* 

4^^2 

k^dk cos 9 d cos 6 

4^^2 



'MaMaM{Et, Ea) - M,MtM{Ea, E,)] , 

'MaMaM{Eb, Ea) + MbMbM{Ea, Eb)] , 
'MaMaM{Eb, Ea) + MbMbM{Ea, Eb)] , 



(41) 



where we have used Eqs. \22\ , (24), and (40), and the isospin indices are chosen to be a = 6 for simplicity. As 
mentioned previously, the longitudinal and temporal components become different to each other, according to the 
breakdown of Lorentz invariance in the temporal direction at finite T. 
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VI. CHARGE SEPARATION, INDUCED MAGNETIC FIELD, AND CENTRALITY IN HIC 



In this Section, we want to discuss ChS defined in Eq. ([3]) using the resuhs given in the previous Section and the 
previous work [15l|T6]. Considering the fact that \{Ji_) Ba,57{J\\}Bo,s\ (x S ^ 1 and (5 <C 1 shown in Ref. [T^, we can 
ignore the term {J±)%^ g rather safely in Eq. (jsj). Moreover, as will be shown shortly, lim|Q|_yo I^^nii^22 c>c X{i,2).Bo.s 
becomes almost zero due to the negative sign between the terms inside the square bracket in the last line of Eq. (41 1. 
Hence, we can drop the term xJ-,-Bo,i5i resulting in the following expression for ChA: 



((cOs(A(/)a + A(/)b))) CX {J3)l„,S + X3,Bo,5- 



Referring to the previous works [151 [TB], {•h)%g & reads: 



k^dk 



2 MaMaM{Ea) 



(42) 



(43) 



,2 _ r2 /12 

''X — " "^0 

which stands for fluctuations in the external-static EM field as in Eq. il7h. We also included the three-dimensional 



where we have used the chiral-chemical potential fj,^ = S^Aq |TS]. Note that we employed Aq, being instead of 1^4, 



Matsubara sum is written as follows: 

M{Ea 



volume V3, resulting in that Eq. (43) becomes dimensionless [10]. The relevant function in Eq. (43 1 coming from the 

1 (l-e-^»/^) 



1 



1 



Ea (1 + e--Ea/T) 2T (1 + e-EjTy 



(44) 



As mentioned, VCC in Eq. (41) in the limit of \Q\ — ^ can be related with the background contribution or the 
susceptibility: 



N 



f 



ImnaadQI) 



(45) 



Here we took an average over the flavor. Again, Eq. (45) is dimensionless. More explicitly, Eq. (451 can be evaluated 
further and resulted in: 



-2n^BoV3NcT lim 



k^dk 



k'^dk cos 9d cos 6 
47r2 

2 MaMaM{Ea,Ea), 



[MaMaMiEb, Ea) + MbMbM{Ea, Eb)] 



where 



MiEa,Ea) 



32T2£;5 



2n 



6Tnanh 



Ea 

2T 



EaSech 



Ea 

2T 



3T + Eata.nh 



Ea 

2T 



(46) 



(47) 



In deriving Eq. (46), we have assumed the following: As seen in Eq. (41), the longitudinal component of connected 



VCC contains a term \Q\5. Considering the relation between 5 and the chiral-chemical potential as above, we can 
write \Q\S = ^^\Q\/Aq. Since the fluctuation Aq must be far smaller than the model scale ^ 600 MeV and can be 
considered nearly zero, we assume the following relation in the limit of (IQj, Aq) — )• 0: 



1. 



(48) 



Due to this simple assumption, one is lead to \Q\5 ^ As a result of all the ingredients disc ussed so far, the relevant 
terms for ChS, {Js}^^ g and X3,Bo,s are expressed as functions of Bq, ^^'i ^ ™ ^1^- (^3) and (46). 

Now, we are in a position to have a geometric consideration of HIC to derive the strength of t tie magnetic field 
as a function of centrality (oriVpait), so as to estimate ChS. A schematic figure for a simple geometry of HIC on the 
transverse plane is depicted in Fig. [3j Note that, since we have set the magnetic field along the z axis, the transverse 
plane is so defined by the x-z plane, being different from the usual definition exploited in HIC, in which the x-y plane 
is chosen for it. As shown in the figure, the radius of the nucleus is assigned as Vnuci, while the normal distance 
between the center of the nucleus and the imaginary center line of the collision as 6/2, in which b denotes the impact 
parameter. For later convenience, we define the region of the overlap of the two nuclei as overlap region. In Ref. [S], 
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the magnetic field induced along HIC was given by the Lienard-Wiechert potential at (r — Tq), where Tq indicates the 
position vector of the charged particle on the transverse plane: 

^/ ^ ^ Z aEM sinhy [(r - rp) x y] 

[{r - ro)2 + {t sinhr - y coshY)^^^^ ' 

where Z, uem, and Y denote the electric charge of the nucleus, fine-structure constant, and rapidity, respectively. If 
we assume only the very early stage of the collision, t Ri y Ri 0, the magnetic field reside on the transverse plane with 
a small rapidity Y <^ 1. For instance, the beam rapidity Ybcam for the Au+Au collision is about 5.371 for y/s = 200 
GeV per a nucleon pair in the nucleus, using the following relation 



y = - In 
2 



PO+Py 
PO-Py 



(50) 



where po and Py is the energy and momentum along the beam (y) direction. Then the rapidity Y can locate at 
—Ybcam <Y< ybcam- Considering the above assumption, the magnetic field can be simplified as follows: 

where we have chosen r — z and Vq = {b/2)x, since we focus on the magnetic field induced on the z axis and 
assumed that the distribution of the magnetic field along the x axis is small. It is worth mentioning that there are 
two electric-charge sources for the magnetic field in HIC: those from the participants inside the overlap region and 
spectators outside region [5]. However, at the very early stage of the collision, the contributions from the two charge 
sources can be nearly the same. In other words, sinhy varies slowly. Moreover, similarly, the electric-charge sources 
are averagely locating at a; = ±b/2 for the early stage as in the assumption. Now, we define a useful quantity, which 
is assigned as: 



Abound = V^nucl " (VS)^. (52) 

As understood by seeing Fig. |3j this lengthy quantity bounds the z value of the overlap region. As discussed in 
Ref. [5 , the magnetic field is screened so that ChS develops mostly from the surface of the overlap region. At the 
same time, this screening effect makes ChS for SCC is larger than that of OCC as discussed in Ref. [5,. Hence the 
screening effect can be parameterized by 

Fscr,±± = exp [-2Q!scr I ^;bound - zl] + exp [-2Q;scr Abound + zl] 

FscY,±T = 2 exp [-Ofscrkbound " ^1] exp [-Q!scrkbound + z\] , (53) 

where ascr indicates a parameter in GeV, which relates to the screening length of QGP: agcr oc 1/Aqgp. Considering 
all the ingredients discussed so far, we can write ChS for SCC and OCC in the following expressions [5]: 

(•Abound ^2 

((COS(A0± -f A0±))) cx C±± / ^ ((J3)lo,5 + X3,Bo,5)^scr,±±, 



-Zbo 



dz 

((COS(A0± + A0^))) cx C±^ / — (( J3)|„> + X3,Bo,5) ^;cr,±T- (54) 

Note that the magnetic field and screening effect are integrated over the overlap region along the z axis. Here Cat is a 
coefficient as a function of the positive- and/or negative-charge particles observed in HIC experiments. For instance, 
it reads ^: 

To make a problem easy, we set C±± sa — C^ip for brevity hereafter. The three-dimensional volume appearing inside 
(Js)^^ g and X3,-Bo.<5 in Eqs. (43 1 and (45) can be estimated for the overlap region as follows: 



b_ b 

V3 « TTZbound^^'nucl = TT^^bound ^ ^nucl • (56) 

Note that Va above is a disc with the overlap area Trzbound^/S and height Tnucb which is a Lorentz-contracted radius 
of the nucleus in the y direction, i.e. 7rnuci- The Lorentz-contraction factor 7 becomes about 100 for -y/s = 200 GeV. 
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FIG. 3: Schematic figure for the transverse plane of a hevay-ion coUision. rnud and b/2 indicate the radius of the nucleus and 
the distance between the nucleus center and the z axis. The magnetic field is generated along the z axis. Thus, the case with 
6 = stands for a head-on collision, whereas h — 2rnuci indicates the case without collision. 



Now we want to discuss the domain with Qt ^ for CME produced in HIC. Since the size of the domain does not 
depend on specific centralities for collisions [3 H] , the probability for finding the domain in a certain volume of 
QGP created in HIC can be regarded proportional to the size of the nucleus, which may corresponds to multiplicity, 
inversely. Similarly, the probability for that the magnetic field induced at the domain can be the same. Therefore, 
the volume in Eq. (56) and the magnetic field in Eq. (51 1 are assumed to be modified as 



Bn 



Bn 



(57) 



In this way, we put the effect of multiplicity for the different nuclei collisions to ChS. Later, we will see that these 
modifications play an important role to reproduce the correct strength hierarchy of ChS for different kinds of nuclei 
depending on their size and atomic mass. 

The next step is to define central! ty with iVpait, which is represented by the number of the participants in HIC, as 
a function of b. Since we already have an expression for the induced magnetic field as a function of b as well as z in 
Eq. (51 1, by doing this, we can relate the magnetic field with centrality of the collision. Again, accounting for the 
collision geometry given in Fig. [3j the centrality A^part can be parameterized with b as follows: 



nucl 



1 - 



tan 



2^7 



nucl 



(6/2)2 



2ri 



V^nucl-(V2)^ 



(58) 



where A'nuci is the atomic mass of the projectile nucleus in HIC. It is easily shown that A'part — 2Nnuc\ for the head-on 
collision {b = 0) and A^part — for the case without a collision (6/2 = fnuci)- For instance, the radius of a gold nucleus 
is about Tnuci ~ 7.27 fm with A^au ~ 197 simply using the well-known relation rnuci — 1.25 fm x N^^^^. Thus, 6/2 can 



vary from to about 7.27 fm for the Au- 
centrality defined as 



Au collision. Then, numerical results for ChS will be given as a function of 



Centrality 



_ (2jVNucl - A^part) 

2N 1 



(59) 



VII. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we provide numerical results for the physical quantities discussed in the previous Sections. First, 
we consider connected VCC defined in Eq. (41). In the left panel of Fig. |4] we show the numerical results for ImHaa 
as a function of \Q\ for different T — (0, 50, 100, 150) MeV. For simplicity, we divided it with 6 and Bq, since VCC 
is linear in these two quantities as in Eq. (41). Thus, one can obtain it for different values of S and Bq by a simple 
scafing. 

As in the figure, at T = 0, we obtain a smooth curve, which has a wide bump in the vicinity of \Q\ = (300 ~ 400) 
MeV, which presents a vector meson corresponding to the current in Eq. ( 11 ). Since we have picked up only the terms 
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proportional to 0{Bq) and 0{6) to examine CME in the present work, the bump position does not match with possible 
physical vector-meson mass, such as p(770) for example. As T increases, the magnitude of the curve becomes larger 
and the bump position shifts to a higher value of \Q\. Moreover, the bump becomes more narrowed. For instance, 
we observe a bump at \Q\ « 600 MeV for T = 150 MeV. We note that this enhancing behavior of the magnitude 
depends much on the lowest frequency wq included in as well as M^. This indicates thermal modification of the 
corresponding vector- meson mass. In Ref. |14j . VCC and corresponding vector-meson spectral function were obtained 
using a SU(2) quenched lattice QCD simulation under a strong magnetic field, and their results for the vector-meson 
spectral function indeed showed the same tendency with that shown in the left panel of Fig. [4] On the contrary, the 
bump position in their work locates at about 1 GeV for T = 0, which is higher than ours and may correspond to 
the mass of p meson. Although we did not show the transverse component of VCC, it behaves similar to that of the 
longitudinal ones but in different magnitudes. 

Here is some physical interpretations on the behaviors shown in the figure. We already mentioned that it depends 
much on the Matsubara frequency inside Et and Mf,. In addition, as noted previously and in Ref. ES], the 
instanton contribution gets diminished as T increases, being indicated by decreasing average instanton size and 
increasing average inter-instanton distance shown in the right panel of Fig. [2j i.e. diluting instanton medium. As 
understood in Eq. (22 1, the relevant part of quark propagator for understanding the problem can be simply written 
as follows: 



{k - qY + A/2 fc2 - {2ttT)\Q\ + |0|2 + Af2 ' 



(60) 



where we have approximated k^^ with the lowest Matsubara frequency wo = ttT. In the light-hand side of the above 
equation, it can be easily noticed that the term (27rr)|(5| plays a kind of a thermal mass, at which a bump is produced, 
while the term Af2 gets smaller due to diluting instanton. Moreover the term {2'kT)\Q\ gets larger linearly with respect 
to T, being consistent with the observation. We verified that, if we choose a higher Matsubara frequency, the bump 
position appears at higher \Q\ values, as expected. For instance, we observe a bump at 400 MeV with wq and at 600 
MeV with wi for T = 50 MeV. At the same time, the magnitude of VCC become enhanced for higher T, since the 



denominator of Eq. (60) gets smaller, being consistent with that of the lattice QCD simulation [T3]. The narrowing 
bump with respect to T also can be explained by decreasing instanton contribution, resulting in smaller M in Eq. ( 60 1 . 

In the right panel of Fig. |4j we show VCC in the limit of |(5| — >■ as a function of \/Bq. Note that this quantity 
is proportional to the susceptibility X3,Bo,5 as in Eq. (45 1. In the calculation, we took \Q\ = 1 MeV as a small value 
to see the behaviors of VCC in the limit of \Q\ — ^ for numerical convenience. Otherwise, one can use Eq. (46 1 for 
correct limiting behaviors. As shown in the figure, the longitudinal component turns out to be finite and increases as 
the magnetic field gets stronger as easily understood by Eq. (|41|). In contrast, the transverse one uniformly zero due 



to the cancelation between the two terms inside the square bracket in the right-hand side of the first line of Eq. (41 1, 
since Ma equals Af^ and, subsequently, M and M are the same for g = 0, equivalently |Q| = 0. In other words, the 
transverse component at 0{Bq) must be zero for all the values of Bq at \Q\ — ^ 0. 

Again, this observation is well consistent with that of Ref. [14] , in which VCC is represented by the corresponding 
conductivities a^z and Uxx- However, there appear differences between theirs and the present calculation as T goes 
higher. The lattice result showed almost flat but nonzero values for the longitudinal and transverse components of 
VCC for higher T. We interpret this difference by that there can be thermal enhancement of VCC for 0{Bq) for 
n 7^ 1 which has not been included in the present work. In our calculations, as for the longitudinal component of 
VCC, one finds that it becomes insensitive to Bq or, equivalently, the slope of the curve gets smaller with respect 
to Bq for higher T, being consistent with the previous work [lSl[T2]. This behavior can be understood as follows: 
Since CME is an effect depending on the topological behaviors of the nontrivial QCD vacuum. Hence, if the vacuum 
contribution decreases due to some reasons, CME also does. We already emphasized that the instanton contribution, 
which represents the nontrivial vacuum, gets weaker and diluted as T increases. As a result, CME and corresponding 
VCC under the strong magnetic field becomes weak with respect to T. In this sense, although we do not go beyond 
in the present work, considering the tendency shown in the figure, it is easy to imagine that the slope of the 
longitudinal component of VCC goes to zero for considerably high T, then becomes flat for Bq finally. A consistent 
situation was observed in the lattice simulation, but with a finite magnitude [14| . This finiteness (nonzero), being 
different from ours, again can be explained by the thermal enhancement 

Now, we are in a position to compute ChS from the instanton contribution at low T defined in Eq. ( 54 ) using the 
present and previous [T3] results for connected and disconnected VCC. Then, ChS can be written Explicitly as follows: 



{{cos{A(l)a + A(f>b))) cx2C„fc 



dz 



k^dk 



1 2 



MaMaMiEa) 



-CT 



k^dk 



MaMaMiEa, Ea) } F, 



sci\ab 7 



(61) 



where we have simply ignored the transverse contributions of VCC as already mentioned and defined a coefficient C 
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FIG. 4: (Color online) ImUss/lSBo) as a function of jQj (left) and li-aUss / {5\Q\) as a function of at \Q\ = 1 MeV (right) 
for different T. In the right panel, we also draw Imnii,22/(<5|(5|) for T = in a thick-solid line. 



as a function of z for convenience: 



C{z) 



fi^Bo{z)V3N, 

nucl 



(62) 



Since we have no available experimental information on fj,^ , we set it as /x^ = 1 MeV as a trial, considering that P- and 
CP-violation strengths are small. Moreover, the value of Cab is chosen as —1 and -1-1 for SCC and OCC, respectively, 
due to that ChS can be easily scaled for reahstic numbers of Na,b as in Eq. (61 1 . Note that we use \Cab\ = 1 even for 
different nuclei collisions, since the dependence of multiplicity (^ size of the nuclei) has been already included in the 
magnetic field and volume as in Eq. (57). 

Before going further, we first look at the T-dependence behavior of (Ja)^ and x in Eqs. (43) and (46 1. Here, we 
have dropped the subscripts Bq and S from the quantities for simplicity. In Fig. [Sj we depict them as functions of T 
up to 3 GeV, which is obviously beyond the applicability of the present framework though. In computing them, we 
set Bq = and V3 — 47rr^^/3 as a trial. As shown in the figure, ( J3) shows a decreasing curve with respect to T 
as already observed in the previous work |15| : Diluting instanton for higher T. As for %, we find a small bump in 
the vicinity T = (50 ^ 100) MeV, then it decreases. Obviously, the magnitude of x is much larger than that of ( J3) 
over the region, whereas the difference between them is diminished and alsmost disappears at T ss 300 MeV. Since 
the sum of these two quantities is responsible for ChS as shown in Eq. (61), we also show the sum of them in the 



figure. As expected, their contribution to ChS is dominated by {J3) ■ This observation is consistent with that given 
in Ref. [9]. Hence, taking into account for the previous and present results, the instanton contribution for ChS (or 
CME) can be almost inferred from the longitudinal disconnected VCC, ( Js)^ alone. 

In Fig.[6) we draw the curves for ChS defined in Eq. ( [6l| ) as a function of centrality in Eq. ([59]). We show separately 
SCC and OCC for ChS for the Au+Au and Cu+Cu collisions for T = 200 MeV and Y = 0.5. The curves in the 
left panel are plotted with ascr = 0.01 GeV, whereas the right one with = 0.5 GeV. Although we did not show 
the numerical reuslt for different T, as for lower T (higher T), one can easily expect that the curves will increase 
(decrease) according to the T dependence of the instanton contributions shown Fig. [s] First, we take a look on the 
left panel of the figure. In general, it turns out that the strengths of ChS (or CME) gets increasing from head-on 
to peripheral collisions. The reason for this increasing is that the induced magnetic field by the collision becomes 
stronger as b increases as understood in Eq. ( 51 ). However, as the two nuclei passing by each other without a collision, 
in other words b — >■ 2rnuci) ChS decreases drastically, since the overlap region disappears by Eq. (561, i.e. Zbound = 0. 



As for the Au-|-Au collision, we observe that SCC turns out to be considerably larger than OCC shown in the left 
panel. This behavior can be understood by the screening effect given by Eq. (53 1: If a particle drifting by CME 



moves inward to the overlap region, the probability to find it after hadronization is much less than that for a particle 
outward from the region |5j. In this sense, the screening effect controls the difference between the strengths for SCC 
and OCC. We observe a similar tendency also for the Cu-|-Cu collision as seen in the figure. However, we find a 
distinctive differences from those for the Au-|-Au collision: 1) the magnitudes of the curves for SCC and OCC get 
larger and 2) the difference between SCC and OCC is diminished in comparison to those for the Au-|-Au collision. As 
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FIG. 5: (Color online) { J3) , X3, and { J3) +X3 in Eqs. (43 1 and (46 1 as functions of T. Here we use Bq 



and V3 = 47rrAu/3. 



mentioned previously, the magnetic field and volume have been modified, accounting for that the domain for Qt =/= 
is not depending on centrality, with the probability ^ 1/A^nuci- Due to this modification, there is strong suppressions 
for ChS for the heavier (or larger) nuclei collision. This is the reason for 1). In contrast, the reason for 2) is a combing 
consequence of the modification and screening effect. In other words, the smaller difference in SCC and OCC means 
the weak screening effect. Since the radius for Cu is smaller than that for Au, the probability that a particle via 
CME drifts through the overlap region then becomes hadronized gets higher. This tendency is clearly demonstrated 
by comparing the left (ctscr — 0.05) and right (agcr = 0.1) panels of the figure. We verified that, as we take smaller 
values for anucb the difference between SCC and OCC becomes smaller, and vice versa for the larger values. 

Here we make a brief discussion on the the value of ascr- In Ref. 0, it was proposed that ascr — V-^QGP- 
Considering the screening length can be estimated as Aqgp oc l/{gT), where g stands for the strong coupling, at 
T = 200 MeV, the length is in the order of less than one fermi. If we take agcr = 1/Aqgp with agcr = 0.05 GeV and 
0.1 GeV, the screening lengths become about 4 fm and 2 fm, respectively. Thus, these values are relatively larger. 
However, since there must be more unknown factors in HIC than what have done in the present work, the simple 



parametrization for the screening effect in Eq. (531 does not compete fully with the experiment, and one needs a 
tuning for the parameters. Moreover, the collision geometry was oversimplified here so that the larger anuci may be 
necessary to compensate what have ignored in the present work. 

We note that the simple collision geometry taken into account in the present work is also valid for different sources 
of CME, such as the sphaleron, as long as the CME current is a linear function of the external magnetic field. 
Considering that the sphaleron contribution is also linear in Bq as in Refs. [51 [TU], the behaviors of the curves in 
Fig. [6] can be valid for higher T but in different magnitudes. In this sense, if we compare the present results for ChS 
with to the experimental data from STAR collaboration [7j[8|, the curves shown in the left panel of the figure may 
be consistent qualitatively with the data, showing the correct strength hierarchy for SCC and OCC for the different 
types of collisions. Moreover, the shape of the curves are similar to the data. 



VIII. SUMMARY AND CONCLUSION 



In the present work, we have investigated VCC and ChS, which was suggested as an indication of CME. To this 
end, we employed the instanton-vacuum configuration with the T modification using the Harrington-Shepard caloron. 
All the relevant quantities were presented as functions of 5, which amounts the strength of P- and CP violations. 
The external magnetic field was included by the linear Schwinger method. Using all these ingredients, we obtained 
expressions for connected and disconnected VCC. The imaginary parts of connected VCC relating to a spectral 
function for a vector meson as well as conductivity were computed. Considering a simple collision geometry of HIC 
and using the present results for VCC, we calculated ChS as a functions of centrality. Below, we summarize the 
important observations of the present work briefly: 



• In leading contributions, connected VCC turns out to be linear in Bq as well as 5. 
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FIG. 6: (Color online) Charge separation in Eq. ( |61[ ) as a function of centrality for the Au+Au and Cu+Cu collisions for 
a = 0.05 (left) and 0.1 (right). We choose T = 200 MeV and Y = 0.5. 



• The imaginary part of connected VCC shows a wide bump at |Q| = (300 ^ 400) MeV for T = 0. The bump 
structure is enhanced and moves toward a higher \Q\ value with respect to T, due to the increasing thermal- 
mass-like effect. At the same time, the bump becomes sharpened by decreasing instanton contribution. 

• Taking the limit \Q\ — ^ 0, the transverse component of connected VCC disappears by the cancelation between 
relevant terms, whereas the longitudinal one remains finite. As T increases, VCC becomes insensitive to i?o, 
denoting a decreasing CME due to diluting instanton. 

• Considering a simple collision geometry for HIC, and assuming the magnetic field generated in a very early 
stage of the collision and screening effect, we can estimate ChS for the Au+Au and Cu+Cu collisions. The 
strength of ChS for SCC turns out to be larger than that for OCC, due to the screening effect, in accordance 
with experiments. 

• ChS for the lighter nucleus collision is more enhanced than that for the heavier one, since CME is proportional 
to 1/-/Vnucr turns out that the difference between SCC and OCC gets smaller as the nucleus involved 
becomes lighter. However, this tendency largely depends on the strength of the screening effect. 

• This simple collision geometry estimation on ChS is almost independent on which source of CME plays a role, 
instanton or sphaleron, as far as the induced CME current is linear in /i^So as in Eq. ([2]). Therefore, the present 
estimation for ChS can be compared with actual HIC experiments to a certain extent. 

We note that there were several simplifications and assumptions to be addressed more carefully in the present work. 
For instance, we have only picked up the terms proportional to 0(Bq) and 0{S) for VCC. By construction, we cannot 
describe the confinement-deconfinement transition, resulting in that we are reluctant to go beyond the critical T for 
SxSB, fa Aqcd- Moreover, the collision geometry was oversimplified which leads to insufficient considerations 
on the complicated macroscopic behaviors of the QGP matter. We do not determine the screening parameter from 
the present model, and it is treated as a free parameter. Especially, the time evolution of ChS was totally ignored. 
However, all the results given in the present work may describe the essence of what happens for VCC as well as ChS 
in HIC, and provide information on the relatively low-T aspect of CME. In addition, the simple collision geometry 
will provide an almost model-independent estimation on HIC, although it is still far from full descriptions on the 
experiments. More sophisticated treatment for the collision geometry, parametrization of Bq, screening effect, and 
time evolution are under progress, and appears elsewhere. 
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